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Abstract 

The Majorana nature of massive neutrinos will be crucially probed in the next-generation 
experiments of the neutrinoless double-beta {0iy2/3) decay. The effective mass term of this 
process, may be contaminated by new physics. So how to interpret a discovery or 

null result of the Oz/2/3 decay in the foreseeable future is highly nontrivial. In this paper we 
introduce a novel three-dimensional description of |(m)gg|, which allows us to see its sensitivity 
to the lightest neutrino mass and two Majorana phases in a transparent way. We take a look 
at to what extent the free parameters of can be well constrained provided a signal of 

the Oz/2/3 decay is observed someday. To fully explore lepton number violation, all the six 
effective Majorana mass terms {m)^^ (for n,/3 = e,/i, r) are calculated and their lower bounds 
are illustrated with the two-dimensional contour figures. The effect of possible new physics on 
the Oz/2/3 decay is also discussed in a mo del-independent way. We find that the result of |(m)gg| 
in the normal (or inverted) neutrino mass ordering case modified by the new physics effect may 
somewhat mimic that in the inverted (or normal) mass ordering case in the standard three- 
flavor scheme. Hence a proper interpretation of a discovery or null result of the 0z^2/d decay may 
demand extra information from some other measurements. 
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1 Introduction 


One of the burning questions in nuclear and particle physics is whether massive neutrinos are 
the Majorana fermions [1]. The latter must be associated with the phenomena of lepton number 
violation (LNV), such as the neutrinoless double-beta (Oz/2/3) decays of some even-even nuclei in 
the form of (^4, Z) (A, Z -h 2) -h 2e“ [2]. On the other hand, the Majorana zero modes may have 
profound consequences or applications in solid-state physics [3]. That is why it is fundamentally 
important to verify the existence of elementary Majorana fermions in Nature. The most suitable 
candidate of this kind is expected to be the massive neutrinos [ 3 ] . 

However, the tiny masses of three known neutrinos make it extremely difficult to identify their 
Majorana nature. The most promising experimental way is to search for the Oz/2/3 decays. Thanks 
to the Schechter-Valle theorem [5], a discovery of the Oz/2/3 decay mode will definitely pin down the 
Majorana nature of massive neutrinos no matter whether this LNV process is mediated by other 
new physics (NP) particles or not. The rate of such a decay mode can be expressed as 

r0" = G0"(Q,Z)|M0"|"|{m),,|y (1) 

where is the phase-space factor, denotes the relevant nuclear matrix element (NME), and 
(m)gg stands for the effective Majorana neutrino mass term. In the standard three-flavor scheme, 

Mee = + "^2^62 + (2) 

with m- (for i = 1,2,3) being the neutrino masses and being the matrix elements of the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino mixing matrix [ 6 ]. Given current neutrino 
oscillation data [7], the three neutrinos may have a normal mass ordering (NMO) 7X1^ < m 2 < m^ 
or an inverted mass ordering (IMO) m 3 < < m 2 - In the presence of NP, (m)^^ is likely to 

be contaminated by extra contributions which can be either constructive or destructive. While an 
observation of the Oz/2/3 decay must point to an appreciable value of |(m)gg|, a null experimental 
result does not necessarily mean that massive neutrinos are the Dirac fermions because (m)^^ ^ 0 
is not impossible even though the neutrinos themselves are the Majorana particles IHIE]. 

Hence how to interpret a discovery or null result of the Oz/2/3 decay in the foreseeable future is 
highly nontrivial and deserves special attention uni El IE]. In this work we focus on the sensitivity 
of |(m)gg| to the unknown parameters in the neutrino sector, which include the absolute neutrino 
mass scale, the Majorana CP-violating phases, and even possible NP contributions. Beyond the 
popular Vissani graph m which gives a two-dimensional description of the dependence of \{m)^J 
on the smallest neutrino mass, we introduce a novel three-dimensional description of the sensitivity 
of |(m)gg| to both the smallest neutrino mass and the Majorana phases in the standard three-flavor 
scheme. We single out the Majorana phase which may make \{m)^J sink into a decline in the 
NMO case, and show that a constructive NP contribution is possible to compensate that decline 
and enhance \{m)^J to the level which more or less mimics the case of the IMO. On the other 
hand, the destructive NP contribution is not impossible to suppress \{m)^J to the level which is 
indiscoverable, even though the neutrino mass ordering is inverted or nearly degenerate. Given a 
discovery of the Oz/2/3 decay, the possibility of constraining the unknown parameters is discussed 
in several cases. We also examine the dependence of \{m)^^\ (for (a,/3 = on the absolute 

neutrino mass scale and three CP-violating phases of the PMNS matrix f7, and conclude that some 
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other possible LNV processes have to be measured in order to fully understand an experimental 
outcome of the Oz/2/3 decay and even determine the Majorana phases. 

2 A three-dimensional description of |(m)ee| 

In the standard three-flavor scheme the unitary PMNS matrix U can be parameterized in terms of 
three rotation angles {0i2^ ^135 ^ 23 ) three phase angles (5, p, a) in the following way [7]: 

Uel = C 12 C 13 , U^2 = «12Cl3 > 

Ue3 = Si3 e'"/' , U^3 = C 13 S 23 , (3) 

where c- = cos 9- and s- = sinO-j (for ij = 12,13,23), 6 is referred to as the Dirac phase since 
it measures the strength of CP violation in neutrino oscillations, p and a are referred to as the 
Majorana phases and have nothing to do with neutrino oscillations. The phase convention taken 
in Eq. (3) is intended to forbid 5 to appear in the effective Majorana mass term of the Oz/2/3 decay: 

IMeel = |"^lC?2C?3 + m'2Sl2C?3 + • (4) 

The merit of this phase convention is obvious. In the extreme case of the NMO or IMO (i.e., 
rrii = 0 or = 0), which is allowed by current experimental data, one of the two Majorana phases 
automatically disappears from \ {m)^J. Note, however, that 6 is intrinsically of the Majorana nature 
because it can enter other effective Majorana mass terms (e.g., {m)^^ and {m)^^ [HI)- 

A measurement of the Oz/2/3 decay allows us to determine or constrain \{m)^J. So far the most 
popular way of presenting \{m)^J has been the Vissani graph m- It illustrates the allowed range 
of |(m)gg| against rrii or by inputting the experimental values of 9^2 and 9i^ and allowing p 
and (j to vary in the interval [0°,360°). In the NMO case |(m)gg| may sink into a decline when 
nil range 0.0023 eV — 0.0063 eV [15], implying a significant or complete cancellation 

among the three components of |(m)gg|. In comparison, there is a lower bound \{ni)^^\ > 0.02 eV 
in the IMO case, and it is always larger than the upper bound of \{ni)^^\ in the NMO case when 
the lightest neutrino mass is smaller than about 0.01 eV m- This salient feature enables us to 
confirm or rule out the IMO, if the future Oz/2/3-decay experiments can reach a sensitivity below 
0.02 eV. Nevertheless, the Vissani graph is unable to tell the dependence of \{ni)^J on p and a. For 
example, which Majorana phase is dominantly responsible for the significant decline of |(m)gg| in 
the NMO case? To answer such questions and explore the whole parameter space, let us generalize 
the two-dimensional Vissani graph by introducing a novel three-dimensional description of |(m)gg|. 

Fig. 1 is a three-dimensional illustration of the lower and upper bounds of \ {ni)^^\ in the NMO 
and IMO cases. In our numerical calculations we have input the best-fit values of Am^i, 

9 12 Q^nd 9^2, obtained from a recent global analysis of current neutrino oscillation data [16]. For 
simplicity, the uncertainties of these four parameters are not taken into account because they do 
not change the main features of |(m)gg|. The unknown Majorana phases p and a are allowed to 
vary in the range [0°,360°), and the neutrino mass or is constrained via the Planck data 
(i.e., nil + ^2 + ^3 ^ 0*23 eV at the 95% confidence level ini)- Some comments on Fig. 1 are 
in order. (1) The upper bound of \{ni)^^\ is trivial, because it can be obtained by simply taking 
p = (j = 0°. (2) The lower bound of \{ni)^^\ is nontrivial, because it is a result of the maximal 
cancellation among the three components of |(m)gg| for given values of p, a and nii or m 3 . (3) In 
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Figure 1: Three-dimensional illustration of the lower (blue) and upper (light orange) bounds of 
|(m)ee| as functions of the lightest neutrino mass and two Majorana phases in the NMO or IMO 
case. 


the NMO case it is the phase p that may lead the lower bound of |(m)gg| to a significant decline 
(even down to zero). In comparison, \{m)^^\ is essentially insensitive to a in both the NMO and 
IMO cases. (4) The allowed range of |(m)gg| in the IMO case exhibits a “steady flow” profile, 
which is consistent with the two-dimensional Vissani graph. Its lower bound (^ 0.02 eV) appears 
at p = 180° for a specific value of and arbitrary values of a, but a deadly cancellation among 
the three components of |(m)gg| has no way to happen. (5) When the neutrino mass spectrum is 
nearly degenerate (i.e., rrii m 2 — m^ > 0.05 eV), the results of |(m)gg| in the NMO and IMO 
cases are almost indistinguishable. 

The parameter space for the vanishing of \{m)^^\ in the NMO case is of particular interest, 
because it points to a null result of the Oz/2/3 decay although massive neutrinos are the Majo¬ 
rana particles. However, the “dark well” of \{m)^J\ versus the p-m^ plane in Fig. 1 has a sharp 
champagne-bottle profile at the ground. This characteristic can be understood by figuring out the 
correlation between and p from \{m)^J\ — 0. Namely, 


244,0 224 I 244 24 

^1^12^13 + 2mim2Ci2Si2Ci3 cos P + m 2 Si 2 Cis = ^3513 


( 5 ) 


Given the best-fit values of Am^i, 9^2 and 9^^ [16], Fig. 2 shows the p-m^ correlation which 

corresponds to the contour of the champagne-bottle profile of \{m)^J in Fig. 1. One can see that 
the “dark well” appears when p lies in the range 160° — 200° and varies from 0.0023 eV to 
0.0063 eV for arbitrary values of a. Such a fine structure of cancellation has been missed before. 

As a matter of fact, a three-dimensional description of \{m)^^\ against two free parameters is 
equivalent to a set of two-dimensional contour figures which project the values of \{m)^^\ onto the 
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Figure 2: A correlation between and p as constrained by the vanishing of in the NMO 

case, corresponding to the contour of the champagne-bottle profile of in Fig. 1. 

parameter-space planes, if only its upper or lower bound is considered. In order to clearly present 
the correspondence between the numerical result of \{m)^^\ and that of a given parameter which 
is difficult to be identified in a three-dimensional graph, we show the contour figures for the lower 
bound of |(m)gg| on the p-cr, m^-p (or m^-p) and m-^-a (or m^-a) planes in the NMO (or IMO) case 
in Fig. 3 (or Fig. 4). For the sake of completeness, we calculate the contour figures for the lower 
bounds of all the six effective Majorana mass terms defined as 

Ma/3 = + m2U^2U^2 + "i3^a3^/33 > (6) 

where the subscripts a and /3 run over e, p and r. There are at least two good reasons for considering 
(a) only the Oz/2/3 decay itself cannot offer sufficient information to fix the three unknown 
parameters of |(m)gg|; (b) if a null result of the Oz/2/3 decay is observed, one will have to search for 
some other LNV processes so as to identify the Majorana nature of massive neutrinos. The typical 
LNV processes which are associated with {m)^^ include the p~ conversion in the nuclear 

background, neutrino-antineutrino oscillations, rare LNV decays of B and D mesons, and so on 
m- In Figs. 3 and 4 the contours for the lower bounds of \ are presented by gradient colors 

and their corresponding magnitudes are indicated by the legends. In particular, the purple areas 
stand for the parameter space where significant cancellations (i.e., < 10“^ eV) can take 

place. When the mg-associated term of is not suppressed by Sl3 ~ 2%, its lower bound 

becomes sensitive to the Majorana phase a. Hence a combined analysis of the 01/2(3 decay and 
some other LNV processes will be greatly helpful to determine or constrain both p and a. 

3 Limits of 3 and p from a signal of the Oi/2/3 decay 

In the standard three-flavor scheme we have studied the possible profile (especially the lower bound) 
of |(m)gg| against the unknown mass and phase parameters. Inversely, the unknown parameters 
can be constrained if the 01/2(3 decay is discovered and the magnitude of \{m)^J is determined. A 
good example of this kind is the strong constraint on the parameter space of mi and p in Eq. ^ 
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Figure 3: The lower bounds of |(^)c^^| changing with rrii^ p and a in the NMO case. 
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Figure 4: The lower bounds of |(^)c^^| changing with m 3 , p and a in the IMO case. 
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or Fig. 2 based on the assumption \{m)^J = 0, which is more or less equivalent to a null result 
of the Oz/2/3 decay provided the experimental sensitivity has been good enough. So it makes sense 
to ask the following question: to what extent the unknown parameters can be constrained from a 
signal of the Oz/2/3 decay? 

Let us try to answer this question in an ideal situation with no concern about the experimental 
error bars. The first issue is to derive the correlation between rrii (or m 3 ) and p like that given in 
Eq. ^ by eliminating a. Since Eq. (4) can be viewed as an implicit function p — f{m-^a) for 
given values of 6 ^ 2 ^ ^13 and |(m)gg|, one may eliminate a by substituting it with the solution of 
dplda\^* = 0. In this way we obtain the maximum and minimum of p as functions of mp 


cos p. 


max,mm 


"^lCl2Cl3 + "^2512^13 “ {^ 3^3 ± l("i)eel)' 

2 ^4 


2 ^ 1 ^ 2 ^ 12 '^ 12^13 


( 7 ) 


If |(m)gg| vanishes, then it is straightforward for Eq. (7) to reproduce Eq. (5). The maximum and 
minimum of a as functions of m- can similarly be obtained: 


cos (J, 


max,mm 


+ mlsf 2 cf 3 - (micf 2 Ci 3 ^ \{'fn)ee\y 


( 8 ) 


However, a is actually insensitive to as shown in Eig. 1. Hence the constraint on a must 

be rather loose even if the Oz/2/3 decay is observed. Eor this reason we simply focus on the possible 
constraints on p and rrii (or m 3 ) in the following. 

Of course, the value of |(m)gg| extracted from a measurement of the Oz/2/3 decay via Eq. 0 
must involve a large uncertainty originating from the NME while the phase-space factor 

G^^{Q^Z) can be precisely calculated. Following Ref. m, we introduce a dimensionless factor 
F to parameterize the uncertainty of \{m)^J inheriting from that of the NME: F = 
where and stand respectively for the maximal and minimal values of the NME which 

are consistently calculated in a given framework. It is apparent that F > 1 holds, and F = 1 cannot 
be reached until the NME is accurately determined. Given a value of F, the “true” value of |(m)gg| 
may lie in the range [| (m )ee\/^^ IM gg|A/F] [IE]- In our numerical calculation we take F = 1 and 
F = 2 for illustration. Fig. 5 shows the allowed regions of (or m 3 ) and p for a few typical values 
of |(m)gg|. The effect of F can be seen when comparing between the cases of F = 1 and F = 2. 
Two comments are in order. ( 1 ) If \{m)^J\ is vanishingly small (e.g., = 0.0005 eV), p can be 

constrained in the range [140°, 220°] in the NMO case. If a larger value of \ {m)^J\ is measured (e.g., 
0.005 eV or 0.05 eV), the allowed range of p will saturate the full interval [0, 360°). To fix the value 
of p needs the input of m^^. Hence some additional information about rrii from the cosmological 
observation or from the direct beta-decay experiment will be greatly helpful. (2) The situation in 
the IMO case is quite similar: p can be constrained in a narrow range if |(m)gg| approaches its 
minimal value (i.e., 0.02 eV), but it is allowed to take any value in the range [0,360°) if \{m)^^\ is 
much larger (e.g., 0.05 eV). Here again is some additional information about m 3 required to pin 
down the value of p. 


4 Possible NP contributions to 


When a NP contribution to the Oz/2/3 decay is concerned, the situation can be quite complicated 
because it may compete with the standard effect (i.e., the one from the three light Majorana 
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Figure 5: The regions of the smallest neutrino mass {rrii or m 3 ) and the Majorana phase p as 
constrained by an “observed” value of In the NMO case |(m)gg| = 0.0005 eV, 0.005 eV 

and 0.05 eV are taken, and in the IMO case |(m)gg| = 0.02 eV, 0.05 eV and 0.1 eV are input. The 
NME uncertainty is illustrated by F. 


neutrinos as discussed above) either constructively or destructively. If the NP effect is significant 
enough, the simple relation between F^^ and |(m)gg| in Eq. (1) has to be modified. This will make 
the interpretation of a discovery or null result of the Oz/2/3 decay more uncertain. Here we aim to 
study the issue in a model-independent way. Namely, we parameterize the possible NP contribution 
to |(m)gg| in terms of its modulus and phase relative to the standard contribution, without going 
into details of any specific NP model [innso]. 

An interesting and very likely case is that different contributions can add in a coherent way 
so that their constructive or destructive interference may happen [El EH- If the helicities of two 
electrons emitted in the NP-induced Oz/2/3 channel are identical to those in the standard channel, 
then the overall rate of the Oz/2/3 decay in Eq. (1) can be modified in the following way: 

ro. ^ z ) 

, (9) 

where denotes the NME subject to the NP process, m^p is a particle-physics parameter 

describing the NP contribution, and (m)gg represents the effective Majorana mass term defined as 

{mYee = rriiFi + 7112^2 + 'm^UeS + ( 10 ) 

with m^p = m^pM^p/M^^. Unless M^p is identical with like the case of NP coming from 
the light sterile neutrinos m^p generally differs from one isotope to another. Hence using 
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Figure 6: The lower (blue) and upper (light orange) bounds of |(m)gg| as functions of (or m3) 
and |mp^p| in the NMO (or IMO) case. 




different isotopes to detect the Oz/2/3 decays is helpful for us to learn whether there is NP beyond 
the standard scenario, but their different NMEs may involve different uncertainties. 

To see the interference between the NP term m^p = |mp^p|e^^NP and the standard one (m)^^ in 
|(m)gg|, we plot the lower and upper bounds of |(m)gg| vs m^ (or m 3 ) and |mp^p| in the NMO (or 
IMO) case in Fig. 6 . For given values of rrii (or m 3 ) and |mp^p|, the lower and upper bounds of 
|(m)gg| can be expressed as 

|(^)L|upper ^ l^elP + ^2 |f^e2l^ + ^3 + I^NpI 5 

iMeellower = max {o, 2171^ \U,f - | M',, , 2|mNp| - | M Lpper } (^1) 

for i = 1,2,3. These results can be directly derived with the help of the “coupling-rod” diagram 
of the Oz/2/3 decay in the presence of the NP [9j. By setting mj,^p ^ 0, we simply arrive at the 
results of \{m)^J obtained before in the standard three-ffavor scheme m- Some comments on our 
numerical results are in order. 

( 1 ) The parameter space in the NMO case can be divided into three regions according to the 
profile of the lower bound of (a) the region with rrii < 0.001 eV and |mj^p| < 0.001 eV, 

where the NP contribution is negligibly small and thus |(m)gg| approximates to 



(b) the region with rrii > 0.01 eV and |(m) 
lower bound 


|(m)ee| - l(m) 


> 


m'lCi2Cl3 “ 


gg| being still dominant over where \{m) 

\Jml + sl^cls - sfg ; 


/ I 

ee I 


has a 

(13) 


and (c) the region with Im^v^pj being dominant over |(m)gg|, where the lower bound of \{m)'^J is 
simply the value of |mj,^p|. If |mp^p| is comparable in magnitude with \{m)^J of the IMO case 
in the standard three-ffavor scheme, it will be impossible to distinguish the NMO case with NP 
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from the IMO case without NP by only measuring the Oz/2/3 decay. This observation would make 
sense in the following situation: a signal of the Oz/2/3 decay looking like the IMO case in the 
standard scenario were measured someday, but the IMO itself were in conflict with the “available” 
cosmological constraint on the sum of three neutrino masses. Note also that at the junctions of the 
aforementioned three regions, |(m)gg| can be vanishingly small either because \{m)^J and 
are both very small or because they undergo a deadly cancellation. 

(2) The profile of the lower bound of in the IMO case is structurally simpler, as shown 

in Fig. 6. In the region dominated by |(m)gg|, \{m)'^J just behaves like \{m)^J in the standard 
scenario and has a lower bound: 


(m) 


- I Me 


> 


2 2 

miCi2Ci3 



2 2 
^ 12^13 


(14) 


On the other hand, \{m)'^J will be saturated by when the latter is dominant over \{m)^J. At 

the junction of these two regions, (m)gg and are comparable in magnitude and have a chance 
to cancel each other. This unfortunate possibility would deserve special attention if the IMO were 
verified by the cosmological data but a signal of the Oz/2/3 decay were not observed in an experiment 
sensitive to the interval in the IMO case of the standard scenario. 


5 Summary 

While most of the particle theorists believe that massive neutrinos must be the Majorana fermions, 
an experimental test of this belief is mandatory. Today a number of Oz/2/3-decay experiments are 
underway for this purpose. It is therefore imperative to consider how to interpret a discovery or 
null result of the Oz/2/3 decay beforehand, before this will finally turn into reality. 

In this work we have tried to do so by presenting some new ideas and results which are essentially 
different from those obtained before. First, we have introduced a three-dimensional description of 
the effective Majorana mass term \{m)^J by going beyond the conventional Vissani graph. This 
new description allows us to look into the sensitivity of \{m)^J (especially its lower bound) to the 
lightest neutrino mass and two Majorana phases in a more transparent way. For example, we have 
shown that it is the Majorana phase p ^ tt that may make \{m)^J sink into a decline in the NMO 
case. Second, we have extended our discussion to all the six effective Majorana masses \ {Tn)^^\ (for 
(a,/3 = e^fi^r) which are associated with a number of different LNV processes, and presented a set 
of two-dimensional contour figures for their lower bounds. We stress that such a study makes sense 
because a measurement of the Oz/2/3 decay itself does not allow us to pin down the two Majorana 
phases. Third, we have studied to what extent (or m^) and p can be well constrained provided 
a discovery of the Oz/2/3 decay (i.e., a definite value of |(m)gg|) is made someday. It is found that the 
smaller \{m)^J\ is, the stronger the constraint will be. Finally, the effect of possible NP contributing 
to the Oz/2/3 decay has been discussed in a model-independent way. It is of particular interest to 
find that the NMO (or IMO) case modified by the NP effect may more or less mimic the IMO (or 
NMO) case in the standard three-flavor scheme. In this case a proper interpretation of a discovery 
or null result of the Oz/2/3 decay demands an input of extra information about the absolute neutrino 
mass scale and (or) Majorana phases from some other measurements. 

In any case it is fundamentally important to identify the Majorana nature of massive neutrinos. 
While there is still a long way to go in this connection, we hope that our study may help pave the 
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way for reaching the exciting destination. 
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